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Abstract. We study properties of the solutions of a family of second order integro- 
differential equations, which describe the large scale dynamics of a class of microscopic 
phase segregation models with particle conserving dynamics. We first establish ex- 
istence and uniqueness as well as some properties of the instantonic solutions. Then 
we concentrate on formal asymptotic (sharp interface) limits. We argue that the 
obtained interface evolution laws (a Stefan-like problem and the Mullins-Sekerka 
solidification model) coincide with the ones which can be obtained in the analogous 
limits from the Cahn-Hilliard equation, the fourth order PDE which is the standard 
Q\ \ macroscopic model for phase segregation with one conservation law. 
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1. Introduction 



In part I [19] we rigorously derived a macroscopic equation describing the time 
evolution of the empirical average process, i.e. the local density, for a model of 
interacting particles with one conservation law undergoing phase segregation. Here 
we establish several properties of the solutions of a class of such equations, which 
are expected to describe more general microscopic models (see Section 3 of [19]). 
Our emphasis is on the late stages of phase segregation, when the phenomena 
are dominated by the motion of sharp interfaces separating domains of different 
densities. 

Let us briefly recall in an informal way the results in part I. The particle models 
are dynamic versions of lattice gases with long range Kac potentials. By a long 
range Kac potential, we mean that the interaction energy between two particles, 
say between a particle at x and one at y (x and y are both in is given by 
7 d J(7(x — y)), where J is a smooth compactly supported function (J(r) = J(—r)) 
and 7 is a positive parameter which is sent to zero. The equilibrium states for 
these models have been investigated thoroughly ([21], [26], [30]) and have provided 
great mathematical insight into the static aspects of phase transition phenomena. 
The dynamical version of these models, in which each particle jumps at random 
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times from a site of the lattice 7L d to one of its unoccupied nearest neighbor sites is 
sometimes called local mean Geld Kawasaki dynamics. The jump times are chosen 
according to a probability distribution which depends on the particle configuration 
and is reversible with respect to the equilibrium Gibbs measure. To find a hydro- 
dynamic scaling limit, we scale also the lattice spacing with 7 and the time with 
7~ 2 (this is the so called diffusive limit). We then derive a (deterministic) evolution 
law for the macroscopic density p. 

In [19] we argue, but do not prove, that our results extend to the case in which 
there is a local (i.e. short range) interaction between the particles in addition to 
the Kac potential: the local interaction should however be sufficiently weak, more 
precisely the system with the local interaction aione should be in a high temperature 
regime. In the present paper we shall consider the equations expected in this more 
general case. They are given in terms of a second order integrodifferential equation: 



d t p(r, t) = V 



a(p)V 
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(1.1) 



p is the density, a (a function of p) is the mobility and 

Hp)= f fc(p(r))dr+] [ [ J(r-r')[p(r)- p(r')] 2 drdr' (1.2) 

J T d 4 J JT<i X T d 

in which f c is either convex or it has a symmetric double well structure, with minima 
at p~ and p + . The latter will be the relevant case for us: we will call the densities 
p~ and p + the phases of the system. The dependence on the temparature (of o 
and f c ) is suppressed. In probabilistic terms, (1.1) is a law of large numbers for 
the empirical averages over the particle system. Equation (1.1) is in a particularly 
enlightening form: it is the gradient flux associated with the classical local mean 
field free energy functional (1.2) ([26,30]), with a density dependent mobility a. The 
form of equation (1.1) allows us to connect the concepts of critical temperature, 
phases, stable, unstable and metastable states of the model, with the properties of 
the solutions of the evolution equation. The next step is to understand the evolution 
of the boundaries (interfaces) between regions in which the density is close to p^. 

Formal matched asymptotic expansions, in the so called sharp interface limit, 
of the solution of macroscopic evolution equations (see e.g. [24,29,5,6]) have been 
successfully employed to understand the interface motion in several models. By 
sharp interface we mean the limit in which the phase domains are very large with 
respect to the size of the interfacial region: if we denote by L the typical size of 
the domains, we will look for results in the limit L — > 00. The time will have to be 
properly scaled as well, typically as some integer power of L, according to the type 
of initial condition (see Section 3) . 

The general picture that we obtain for the interface motion is the following: 
choose an initial condition which takes only metastable or stable values over large 
domains(of typical diameter O(L)), while the interfacial regions are layers of thick- 
ness O(l) and let it evolve according to (1.1). There is, first, the equilibration of 
the interface which happens on a fast time scale (t <C L 2 ). Then, on the time scale 
t oc L 2 the evolution of the density in the bulks (that is the interior of the domains) 
is given by a nondegenerate nonlinear diffusion equation with Dirichlet boundary 
conditions on a free boundary, the interface (Stefan problem). Once the density 

2 



in the bulks is relaxed to the density of the phases, the motion of the interface 
on this time scale stops. A slower evolution can then be seen on the time scale 
t oc L 3 and the motion of the interfaces is given by the Mullins-Sekerka model; 
a quasistatic free boundary problem in which the mean curvature of the interface 
plays a fundamental role. 

Throughout the paper we will repeatedly stress the fact that the evolution in the 
sharp interface limits we have just outlined are the same as those obtained from the 
Cahn-Hilliard equation in the corresponding limits. The Cahn-Hilliard equation 
[7] is an equation of the form (1.1) but the free energy functional is chosen of the 
Ginzburg-Landau type 



in which V is a smooth function (say a polynomial) with a double minimum since 
the aim is to study phase segregation phenomena. The results on sharp interface 
motion for the Cahn-Hilliard equation have been obtained formally by R. L. Pego 
[29] and part of his program has been made rigorous recently [1,3]. In [29] only the 
case a = 1 is considered, but the general case is analogous as long as the mobility 
does not vanish in the phases (the minima of V). The differences between the 
two models in this case are mostly in their having different surface tension, which 
appears as a factor in the formulas: the parameters of the models can therefore be 
tuned to give exactly the same sharp interface motion. A possibly deeper difference 
arises in the case of zero temperature (in our case the temperature is built in the 
model while in the Cahn-Hilliard case one has to choose an appropriate V which 
is temperature dependent, see e.g. [6]): the profile of the stationary front solutions 
(instantonic solutions) in the two cases has a substantially different limit. The 
instantonic solution in our case approaches a step function, while in the C-H case 
it approaches a differentiable function. The zero temperature case has been treated 
(formally) in [6] , but here we will limit ourselves to a few observations in this case. 

Finally we remark that in the case of particle dynamics without conservation 
law (Glauber or Kawasaki+Glauber dynamics [32]) there are by now several results 
proving that the interface motion is by mean curvature ([2,12,22,33]). 

The paper is organized as follows: in Section 2 we introduce in an informal 
way the equations and all the relevant ther mo dynamical concepts. In Section 3 we 
introduce the sharp interface limit and state the results of the formal asymptotics: 
the computations are carried out in Section 5. In Section 4 we state and prove 
existence and uniqueness for the equations introduced in Section 2 as well as a 
result of existence and uniqueness of istantonic solutions. 



In this Section we introduce the evolution equations in an informal way, assuming 
that all the functions which appear in the formulas are suitably smooth. Moreover, 
in order to understand better the role of the quantities appearing in the definitions, 
we will restrict ourselves for the moment to the cases in which we will study phase 
segregation phenomena. Many of the assumptions made here are not needed to 
establish properties like existence and uniqueness: these will be stated and proven 
in a more general context in Section 4. 




(1.3) 



2. The nonlocal equation 
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Let T d = R d mod aL d be the torus of diameter a G M + (seen as a metric space, 
equipped with the periodic Euclidean distance). Occasionally we will consider the 
case a = +00, that is T d = M d , but, unless esplicitly stated, a has to be considered 
fixed (and finite). 

The free energy and the mobility. We start by introducing the free energy 
functional 

HP)= I f{ P {r))dr-\ f f J(r-r')p(r)p(r>)drdr' (2.1) 

Jrpd Z J Jj'd xT d 

in which p is a function from T d to [0, 1], / : [0, 1] — > R is convex and we will refer 
to f(p) as the free energy of the reference system corresponding to the interaction 
potential J being zero. J is chosen attractive, 

J > 0, (2.2) 

isotropic, 

J(r) depends only on \r\ (2-3) 

and is compactly supported in a ball of diameter smaller than the diameter of T d . 

Define moreover the mobility cr, a function from [0, 1] taking nonnegative values 
and such that cr(0) = cr(l) = 0. 

Further assumptions on / and a. We will make three further assumptions on 
/ and a: 

1. Symmetry. We assume that both / and a are symmetric with respect to 1/2, 

1. e. 

/((I + m)/2) = /((l - m)/2), a((l + m)/2) = <r((l - m)/2) (2.4) 

for all m G [—1,1]. This is a consequence of the particle-hole symmetry in the 
microscopic dynamics. 

2. The reference system is a nondegenerate diffusion. We will require that there 
exists c G [1, 00) such that 

- c <D(p)^a(p)f"(p)< c (2.5) 

for all p G (0, 1). We will assume that the limits of D{p) when p approaches and 
1 exist and these will extend the definition of D(-) to [0, 1]. 

3. There are at most two phases. If we set 

Up) = - J -^(p- 1 ^) +f{ P ) (2.6) 

where J(0) = J J(r)dr, the free energy (2.1) becomes, up to a irrelevant additive 
constant, the functional defined in (1.2). If p is constant the second term in (1.2) 
vanishes, so that f c is the free energy density of a homogeneous profile (or con- 
strained equilibrium free energy, see [26, 30]). For (small enough) values of J(0), f c 
can be convex: in this case there is no phase segregation. We are then focusing on 

4 



the nonconvex case and we assume that f c has a double well structure: precisely, 
there exists p+, p~ G (0, 1), 2p+ —1 = 1 — 2p~ , such that 

/c(p)<0, ifpe(p",p+) 
£'(p)>0, if P e(0,p-)u(p+,l). 

Once / c is non convex, assumption (2.7) is, for example, consequence of the stronger 
assumption 

f"(p) is strictly increasing for p > 1/2 (2-8) 
which holds if the reference system is ferromagnetic (see Section 1 of [25]). 
The evolution problem. We consider the equation 



d tP = V 



(2.9) 



v<w. 

with initial condition p(-,0) = p(-). An alternative way to write (2.9) is 

d t p(r, t) = V [D(p(r, f))Vp(r, t) - a(p(r, f))(V( J(-) * p(-, f))(r)] (2.10) 

where * denotes the convolution in T d . Global existence and uniqueness for a 
weak version of (2.10) is given in Theorem 4.1; see Remark 4.1 for the regularity 
properties. 

Stable, unstable and metastable states. By assumption (2.7) the interval 
(0, 1) can be partitioned into three sets: 

(1) the unstable region 

in which f'J < and p^ are defined by (Pm) = 0; 

(2) the metastable region 

[P~,Pm)U(p+,P + ] (2.12) 

where p^ are defined by f' c (p^) = 0; 

(3) the stable region 

[0,p-)U(p+l] (2.13) 

in which f" is bounded away from zero (and it is +oo at and 1). The two minima 
p ± are called phases. 

All the definitions and assuptions we have made are consequences of the statis- 
tical mechanical origin of the model and they are explained at length in [18] and 
[19]. In particular the fact that the density p is between and 1 follows from the 
original particle model which cannot have more than one particle at a lattice site. 
Another consequence of this is that cr(0) = cr(l) = 0. The free energy / comes from 
a particle model which is above its critical temperature and it is then a smooth 
function in (0, 1) [31]. The control on the regularity of the mobility a is harder, but 
Lipschitz continuity (see (4.8)) can be shown to hold in a rather general context 
[32]. The assumption (2.5) is natural, since D is the diffusion coefficient of a high 
temperature system (see [19] and [32]). We stress however that there is no funda- 
mental reason to assume attractivity and isotropy of J(r) to derive (2.9). Even in 
a non attractive case interesting phase segregation phenomena might occur: this is 
not considered in this paper. Moreover, in the anisotropic case one can do a formal 
asymptotic analysis as in Section 5 and get more general results than the ones we 
present; see (4.9) for a condition which replaces the isotropy. 
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A particular case. For the particle system treated in [19] with complete proofs 
we have free energy density 



f(p) = ^(plogp + (1 - p) log(l - p)) (2.14) 

and mobility 

*°(p) = Pp(l-p) (2.15) 

where j3 is a positive parameter, corresponding to the inverse of the temperature. 
In this case, for (3J(0) > 1/4, f"(l/2) < 0, while for /3J(0) < 1/4 we have that 
f"(p) > for all p G (0, 1). Notice that in this example D = 1. 

3. The sharp interface limit 

We want to study the motion of boundaries between regions in which the density 
is either stable or metastable. Our analysis is based on formal matched asymptotic 
expansions. 

Space— time rescaling. in what follows we shall consider the evolution (2.9) with 
r in e~ 1 T d with e > and we will scale time as e~ q , with q = 0,2,3 (only a 
short remark will be made about q = 4), according to the case. We will make 
the substitution t = e~ q r q . At this stage, it is more convenient to reformulate the 
problem in terms of the rescaled densities 

p £ (r,r 9 ) = / 9(e- 1 r,e-%) (3.1) 

in which the spatial variable r ranges in T d and r q is a positive number. The 
evolution equation (2.9) yields 

e q ~ 2 d Tq f(r, r q ) = V [a(p £ (r, r q )) V/* e (r, r q )] (3.2) 

in which 

^(p £ )(r,r) = /V(r,T))- / J e (r - r')p*(r' ,r)dr' (3.3) 

JT d 

and J e (r) = e~ d J(e~ 1 r). Recall that J is taken to be supported in a fixed ball, 
independent of e. 

Before starting the analysis of the time evolution of the system, we need to 
introduce the concept of instantonic solution. 

The one dimensional problem and its instantonic solution. We replace T d 
with M. d in (2.9) and let the initial condition p°(r) = wo(^i) depend only on the 
1— coordinate. Then (2.9) reduces to the one dimensional problem (r\ = z) 

d t u(z, t) = d z [a(u(z, t)) (d z f(u(z, t))-d z ((J* u)(z, t))] (3.4) 

in which 

J(z)= J(z,r 2 ,... ,r d )dr 2 ...dr d (3.5) 

The initial condition is u(-,0) = uq(-). Clearly u(z,t) = C G [0,1] is a stationary 
solution. It is easy to see that, if C is in the unstable region, this solution is unstable. 
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In Section 4, Proposition 4.2, we show that (3.4) admits a nonconstant stationary 

solution which is unique up to translation and reflections in a proper class (see 

Proposition 4.2 for the precise statement). U solves the functional equation (see 
the proof of Proposition 4.2) 

U(z) = (fT 1 ((j*(u- 1 -yz)^ (3.6) 

in which (/') _1 is the inverse function of /', which is invertible since we assumed 
/" > 0. From (3.6) it follows directly that 

lim U{z) = p ± (3.7) 

2— >±CSO 

hence the asymptotic values for U, which we will call the instantonic solution, are 
exactly the phases. We observe that p(r, t) = U(zq + u • r) is a stationary solution 
of (2.9) in the case T d = R d for any choice of z e R and v e M d , \u\ = 1. This 
solution represents a plane stationary front. 

The sharp interface initial condition. The type of initial density profiles p e Q 
that we wish to consider take values which are either metastable or stable on the 
whole space T d , with the exception of the points which are at distance 0(e) from a 
smooth (hyper) surface To C T d , which will be called the interface. We will denote 
by u(r) one of the two unit vectors normal to the interface at the point r: the 
(conventional) choice of the orientation will be explained below. We suppose that 
there exists a function p , defined in T d and smooth outside r , such that if r is 
not on the interface, then 

lim p e (r) = p (r) (3.8) 

and the corresponding limit holds also for the derivatives of Pq. Moreover if 
4>(r, r ) = 0(e) (4> is the signed distance), then 

pl(r) = U(e- 1 (t>(r,r )) + O(6) (3.9) 

In particular the limit of po(r) as r approaches To from the interior of a cluster is 
either p + or p~ and po is discontinuous at To- The sign of <p(r, To) is chosen to be 
positive (respectively negative) if po(r) > 1/2 (respectively po(r) < 1/2) and v is 
chosen to point in the direction in which <p increases, i.e. v = V0 on the interface. 

It turns out that the relevant time scale for the evolution of To is t oc e~ 2 , that 
is q = 2 in (3.2). 

A remark about the equilibration of the interface. The sharp interface 
initial condition that we introduced is such that the interface is locally stationary, 
since it is locally very close to the plane stationary front. It is however very natural 
to consider initial condition which satisfy (3.8), but not necessarily (3.9). In partic- 
ular, the limiting valuess of po on the interface may not coincide with the phases. In 
this case, the first observation is that d t p e (r, t) = 0(e 2 ) for t = and r away from 
the interface. However if r is close to the interface then dtp £ (r, t) = O(l) at t = 0. 
We then expect the transition layer around the interface to evolve on the time scale 
O(l) and eventually (asymptotically on the time scale 0(1)) to approach a profile 
which, to leading order, is given by U(e~ 1 (p(r 1 To))- This is a strong ansatz on the 
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behavior of the solutions of (2.9). Due to the conservation of mass, an immediate 
problem arises: the two asymptotic values of U, i.e. the two phases, will not in 
general match with the solution away from the interface. This happens whenever 
the values of po near To do not coincide with the phases. In this case we have 
to deal with a boundary layer problem. This early stage of the evolution can be 
studied, at least on a formal level, and it gives rise to a self similar Stefan problem: 
the formal expansion is absolutely identical to that in Section 3 of [29]). We will 
not focus on this stage of the evolution and we will only remark that the boundary 
layer should be of thickness 0(ey/i), so that it disappears if t = 0(e -2 ). We will 
focus on the later stage by assuming that this stabilization has taken place. 

Sharp interface limit I: Stefan problem (q = 2). If we assume that the 
interface is stable (see [1] for some rigorous statements on stability of interfaces 
for the Cahn-Hilliard model) the structure of the initial condition is essentially 
preserved and hence p e (r, T2) (together with its derivatives) converges to p(r,T2) 
away from a hypersurface Y T2 while the boundary values of p(r, T2) at V coincide 
with the phases. Moreover we assume that p is smooth outside T T2 and that T T2 is 
itself smooth and (d — l)-dimensional. In general this will be true only for T2 in a 
finite interval and when a singularity appears one should introduce a weak version 
of the evolution law. We have the natural domain decomposition 

T d = r T2 uo+ uq- (3.10) 

and the sets appearing on the right-hand side of (3.10) are disjoint. 

Statement 1. In the setting specified above and at the level of formal asymptotics, 
P e (r, T2) converges to p(r, T2), which solves the following nonlinear Stefan problem 

d T2 p = V-(a(p)Vp (p)) forren± 

linv^ (r) ± p(r', r 2 ) = p ± for r e T T2 (3.11) 
p(r,0) = p o (r) for all re T d 

in which po(p) = f' c (p) and by limit as r' goes to (r) 1 * 1 we mean the limit A — > 
with r' = r + \v . The motion of the interface is generated by the velocity field 
V2(r)v(r), where r e T T2 and 

V ^= {2p f-l) fr'^-fr) ( 3 - 12 ) 

in which [v-'Vpo]_(r) = lim r /_^( r )+ u(r) ■ Vpo(r') — lim r /^( r )- u(r) ■ Vpo(r'), i.e. the 
jump in the normal gradient of po across the interface and a ± = cr(p+) = o~(p~). 

Remark 1: in the case Tq = 0, that is the density takes only stable or metastable 
values and there is no interface, this result has been made rigorous, even starting 
directly from the particle system (see references in [19]). 

Remark 2: The Stefan problem we just introduced is expected to stabilize as 
T2 — > 00 and it should lead to homogeneous density in the (surviving clusters). This 
fact is however very nontrivial from the mathematical viewpoint and the problem 
(3.11), (3.12) may become ill-posed for some finite T2. It is however easy to verify 
that if the density in the clusters is that of one of the homegeneous phases, then 
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V 2 = and so there is no macroscopic motion of the interfaces on the time scale 
q = 2, but there may be a motion on a longer time scale, i.e. if q is larger. 

The sharp interface II: the Mullins— Sekerka symmetric solidification 
model (q = 3) . Let us assume that the initial condition satisfies the extra condition 
Po( r ) = P ± + 0( e ) f° r every r away from the interface and that, like for the case 
q = 2, the density profile at the interface is given by the rescaled instantonic 
solution in the direction perpendicular to the interface. It turns out that, for such 
an initial condition, the relevant time scale is q = 3. Under the assumption of the 
stability of the interface (see the discussion for the q = 2 case) , we can describe the 
time evolution of the density in terms of the evolution of the interface, since the 
density in the cluster is just the equilibrium one. The problem is then reduced to 
find the geometrical motion of the interface r T3 . The analog of the natural domain 
decomposition (3.10) clearly holds also in this case. In the Appendix it is shown 
that for this model the surface tension (i.e. the free energy density associated with 
the instantonic solution U, see [4]) is 

S=-[ {z - z)d z U(z)U(z')J(z- z')dzdz' (3.13) 
2 JrxR 

Let us furthermore denote by K(r), r on the interface, the sum of the principal 
curvatures of the interface at r, that is d — 1 times the mean curvature at r. 

Statement 2. In the setting explained above and at the level of formal asymptotics, 
the motion of the interface Y Tz is generated by the velocity field V^(r)v(r), for r 
belonging to the interface Y T3 . The speed is given by the following expression 

Vs(r) = ^- T [u-Vfi 1 ] + (r) (3.14) 

where [v ■ V p,\]_(r) denotes the jump of the component normal to the interface of 
the function which is continuous on T d , at least twice differentiable on and 
solves the static problem 



A//i(r) = /orrefi+un- 
^(r) = -SK(r)/(2p + - 1) for r £ T T3 [ ' } 

Remark 1: The global existence for the Mullins-Sekerka evolution law presents 
severe mathematical difficulties. We refer to [8] for these aspects of the problem. 

Remark 2: The function pi, which appears in (3.14) and (3.15), may be thought 
of just as an auxiliary tool for computing the interface velocity. It does however 
have a direct physical meaning: it is the 0(e) correction to the chemical potential 
away from the interface, i.e. (i = nip^) + epi + o(e) in the interior of the clusters. 

Some remarks on the limit of vanishing temperature. Here we restrict 
ourselves to the model defined by (2.14) and (2.15). To stress the dependence on 
the parameter f3 we add a subscript (3 to U, a ± , p ± and S. In this natural 
question is what happens to the velocity V3 in the Mullins-Sekerka model as (3 
goes to infinity. In this limit the instantonic solution Up approaches the Heaviside 
function (easily shown using formula (3.6)), i.e. the characteristic function of the 
positive semiaxis. Therefore, using (3.13), we easily see that 



1 [°° ~ 
lim Sp = Soo = - / rJ(r)dr 



and is strictly positive and finite. By (3.15) the gradient of p,\ is uniformly 
bounded as long as the interface is smooth. So what really matters to our analysis 
is the behavior, as j3 becomes larger and larger, of the factor /3cr^/(2p^ — 1) in 
(3.14). While the denominator tends to one, the mobility in the phase ap vanishes 
exponentially fast. Our conclusion is then 

lim V3 = 

and the limit is approached exponentially fast. So the Mullins-Sekerka motion is 
rapidly frozen when the temperature is lowered. A motion on the time scale q = 4 
is expected to arise (surface diffusion), but this will not be considered here. 

4. Existence, uniqueness and regularity 

We follow closely the notation of [10], to which we refer for all the basic notions 
and results. In the first part of this Section we work in T d (unit torus), so that 
L 2 = L 2 (T d ), H 1 = H 1 ^), .... r > is fixed and our basic (Hilbert) space is 



W = W(0,t; H 1 ,!!- 1 ) = : w G L 2 (0, r; H 1 ), ^ G L 2 (0, r; H' 1 ) j (4.1 



) 



where H~ l = (H 1 )' , the dual of H 1 with respect to the L 2 scalar product (•,•)• 
The norm in H 1 will be denoted by || • ||, the norm in L 2 by | • I2 and the norm in 
H~ x by || • II*. In particular we have 



in which v G H~ x and v is the Fourier transform of v. The space W is equipped 
with the (Hilbert) norm 

/ \ i/ 2 

\H\\V = (\\w\\h(0,T;Hi) + \\ W 'Wh(0,r-H-^) 

T 1/2 

= [\Ht)\\ 2 + \\w'(t)\\l]d?j (4.3) 

in which w' = dw/dt. We will often work with the following convex subset of W 
W 1 = {w e W : < w(t, x) < 1 for almost all (t, x) G [0, r] x T d }. (4.4) 
We say that p G W\ is a weak solution of (2.9) if for all u G H 1 

^ («, p) + (V«, D{p) Vp - a(p) (V J * p)) = (4.5) 

in the sense of £>'((0,t)) (X>((0, r)) is the space of C°° compactly supported real 
functions over (0, r)) and 

p(0) = p. (4.6) 
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We recall that (4.6) makes sense because W can be continuously imbedded in 
C°(0, r; L 2 ), Th.l, Ch.XVIII of [10], and that d {u, p) /dt = (u, p'). 

The theorem we are going to state below holds under the assumptions made 
above and under the following assumptions on D, a and J. 

DeC°([0,l]) (4.7) 

(u E C k ~ , k E Z + , means that the derivatives of u of order (k — l) th are Lipschitz) 
and that 1/c < D < c for some c > 1. Moreover 

a E C 1_ ([0, 1]) a{p) > for all p E [0,1], <x(0) = <r(l) = 0. (4.8) 

Finally 

j e C 2^ T d^ = j(_ r ) for aU r e T d_ ^ 

Theorem 4.1. There exists a unique solution to the problem (4.5), (4.6). 
Proo/. 

Existence. We will make use of the following result. 
Proposition 4.1. For every g G £ 2 (0, r; L 2 ) £/ie problem 

u' = V{D(u)Vu - <r(u)(VJ * g)} 



(4.10) 



u(0) = p 

/ias a unique weak solution (in the same sense as (4.5) and (4.6),). 



Proof of Proposition J±.l: (4.10) is a parabolic problem for which existence in 
W can be established via fixed point theorems, see e.g. [27], IV. 4, (extending 
the definitions of a and D outside [0, 1] in a arbitrary (nice) way). The comparison 
principle and the fact that <r(0) = cr(l) = 0, so that u = 1 and u = are solutions of 
(4.10) with proper initial condition, immediately allows to obtain a solution in W\. 
Uniqueness is also standard and can be, for example, established along the same 
line that we are going to use below to prove uniqueness for the original problem 
(4.5), (4.6). □ 

Proposition 4.1 defines a map X : L 2 (0, r; L 2 ) — >■ L 2 (0, r; I?) 

X(g) = u (4.11) 

and we notice that, since < u < 1 almost everywhere, ||X(g) ||| 2 ( T . L ^ < r. 
Directly from (4.10) one has that 

(u(t),u'(t))+ [ \Vu(t)\ 2 D(u(t))= [ o-{u{t)){VJ*g{t))Vu{t) 

J T d J T d 

so that by using the Cauchy-Schwarz and Young's inequalities, (2.5) and the fact 
that a is bounded, one obtains that there is a finite constant c\ such that 

^l«(*)H + ^|V«(t)|5<y 1^)11 (4-12) 
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in which c is the constant in (2.5). From (4.12), by integrating in time, we extract 
the following bound 



\u(t)\\ 2 dt < c 2 



L2( ,r;L2) 



+ 1 



(4.13) 



for some finite constant c 2 (independent of the initial condition). Moreover 



I it' 1 1* = sup I / [ff(u)(VJ * fiOViti — D(u)VuVu] 



which by the Cauchy-Schwarz inequality implies that there is a cs such that 



I«'(f)ll2d*<<* 



I i|2 .-I 
\9\\l 2 (0,t;L 2 ) + 1 



(4.14) 



so that by (4.13) and (4.14) we have that there exists C4 such that 



\u\\w < C4 



1 11 2 , -1 

\9\\l 2 {o,t-l' 2 ) + 1 



(4.15) 



which implies that X is compact, since W is compactly imbedded in L 2 (0, r;L 2 ), 
see Proposition 4.2, Section IV of [27]. X is also continuous since given a sequence 
{g n } which converges in L 2 (0,r;L 2 ) to g, the sequence {X(g n )} = {u n } is rel- 
atively compact in L 2 (0, r; L 2 ) and weakly relatively compact in W (by (4.15)). 
This in particular implies that D{u n )Vu n converges weakly along subsequences in 
L 2 (0, r; L 2 ) to D(u)Wu, since D(-) is C° and bounded. Hence any limit point u of 
{u n } (weakly in W and strongly in L 2 (0, r;L 2 )) solves (4.10), which has a unique 
solution, and the continuity of X follows. 

Summing up, we have obtained that the continuous map X takes the closed 
convex set 

{ueL 2 (0,r;L 2 ) : |M|£ 2(0)T;L2 ) < r} 

into itself and the image of this set under X is relatively compact (in particular 
X(u) G W\ and ||w||w ^ c 4[l + r1 ^ 2 ])- By the Schauder fixed point Theorem 
([20], Corollary 11.2) X has a fixed point p £ W\ which satisfies (4.5), (4.6). This 
completes the proof of the existence. 

Uniqueness. 

Call pi and p 2 two solutions of (4.5) such that pi(0) = p 2 (0). Set X(p) = 
f£D(p')dp', so that A e C 1 and X'(p) > 1/c > 0. We have 



1 d 

2dt 



Pi(t) ~ P2(t)\\l = - (V( Pl (t) - P2 {t)\ (1 - A)- 1 V(A( Pl (t)) - A( P2 (t)))) + 



(V(pi(t) - /9 2 (t)), (1 - A)" 1 V(a(p 1 )VJ * Pl (t) - a(p 2 )VJ * p 2 (t))) . 
We then have 
1 d 



2dt 



\ P l(t) ~ P2(t)\\: = - (Pl(t) - P2(t), (X(pi(t)) - X(p 2 (t)))) 

+ ( Pl (f) - p 2 (t), (1 - A)-\X( Pl (t)) - X(p 2 (t)))) 
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+ (V(pi(t) - p 2 (t)), (1 - A)" 1 V(a(p 1 )VJ * Pl (t) - a(p 2 )VJ * p 2 {t))) . 

Since A' > 1/c > and by the Young and Cauchy-Schwarz inequalities, we obtain 
that for any €1,62 > 

^IM*) " P2<t)\\l < ~\ |Pl(*) - ^(^ + |||A(P1) " A( P2 )||2 + ^2 IIPI - 



1_„„, ,,,0 e 2 

Since for all v G L 2 , ||v|| + < \v\ 2 , for all it G if 1 , ||Vw||* < |tt j 2, and 



+ 2^l|V(pi - p 2 )\\l + 4HVk(pi)VJ * Pl (t) - a{p 2 )VJ * p 2 (t)]\U- (4-16) 



||V[(r(pi)VJ*pi(t) -(7(p 2 )VJ*p 2 (t)]|L < \Wpi)VJ * pi(t) - a( P2 )VJ * p 2 (t)\ 2 

< \ W(pi) - a(p 2 )]\VJ * pi\\ 2 + \a(p 2 )\VJ * (p 2 - Pi)\\ 2 < c 5 \pi - p 2 \ 2 

where C5 depends only on the Lipschitz constant of a and on |VJ| 2 . From (4.16), 
choosing e\ and e 2 small, we obtain that there is cq such that 

\^ t \\pi(t)-p 2 (t)\\l<C6\\pi(t)-p 2 (t)\\l 

which, since ||pi(0) — p 2 (0)||* = 0, implies that ||pi(t) — p 2 (t)\\* = for all t G [0,r] 
and, since pi,p 2 G W, \\pi(t) — p 2 (t)\\ = which implies the result. □ 

Remark 4.1: We observe that since p G C°(0, r; L 2 ), VJ * p G C°(0, r; C 1 ) (since 
J G C 2 ). Once we have a solution of (4.5), (4.6), we have also a solution of (4.10) 
for g = J * p. It is then a standard result (see e.g. [23]) that, if p G C 2 and 
if D and a are C 1 , (4.10) has a classical solution, which is unique in the class of 
functions considered in Theorem 4.1 and Proposition 4.1. Therefore under the same 
assuptions one has a classical solution for our original problem (4.5), (4.6). One 
can obtain further regularity properties by making further assumptions on D, a 
and J. 

We now establish the existence and uniqueness (in a proper sense) of the is- 
tantonic solution and we give some of its properties. Sometimes in what follows 
we keep the notation with partial derivatives also for functions of one variable, for 
uniformity. 

Consider the one dimensional setting with a = 00, i.e. we are working on M. 
Assume (2.2), (2.4), (2.8) , that f c is non convex and all the assumptions made 
for Theorem 1. For consistency with the sharp asymptotic computations, below we 
use J for J in the one dimensional setting (recall (3.5)). Take J to be compactly 
supported. 

Proposition 4.2. Under the assumptions stated above, there exists a function U G 
C^R) such that for all v eC^< 



/ 

Jr 



d z v{z) D(U)d z U(z) - a(U(z))d z (J *U)(z) dz = 0. (4.17) 



U satisfies: 

(1) U - 1/2 is odd; 
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(2) U'{z) > for all z G R; 

(3) There exists m > such that 



lim e m|z| \U(z) -p ± \ = 0. (4.18) 

Moreover if a non constant non decreasing function u G H} oc (M), u(z) G (0, 1) for 
all satisfies (4.17) with U = u for all v G Cg°(R), then there exists I6l 

smc/i. £/ia£ it(z) = U(z — z). 

Proof: take it as in the statement of the Proposition above. Then there exists 
C 6R such that 

D{u)d z u{z) ~ o-(u(z))(J' * u){z) = C (4.19) 

for almost all z. Since the second term in the left-hand side vanishes as z — > oo 
and D > 1/c, then C > and d z u(z) > cC/2 for 2 sufficiently large, so that C = 0. 
Since it(z) G (0, 1) for all z, \f"\ (recall (2.5)) is bounded on every compact subset 
of R and a(u) > 0. Hence we have 

d z [f'(u(z))-(J*u)(z)] =0 (4.20) 

which, using the fact that /" > 0, is equivalent to 

u(z) = {f')- 1 ((J * (« - l/2))(z) + fc) (4.21) 

for some h G R. In [16] and [14] it is shown that, in the class of functions we are 
considering and for the particular case f'{p) = log(p/l — p), so that (f')~ 1 (x) = 
(l/2)[tanh(x/2) + 1] there is a solution to (4.21) in the case h = and [14] it is 
unique in a class larger than the one we consider. Their approach can be repeated 
in our context, but the general result in [9], in particular Remark 5.2 (ii), keeping 
in mind that we have the symmetry condition (2.7), covers the case considered 
here, establishing thus existence and uniqueness for h = 0: property (1) follows 
from (2.7) as well. Property (2) is established by the argument on page 706 of [14] 
and property (3) is established by repeating the proof of Proposition 2.2 in [13] 
for the general / considered here. We have thus established the existence of the 
function U claimed in the Proposition and that U solves (3.6). We remark that 
in [11,14,13], J(0) = 1 and J is multiplied by f3 > 0, while in our case the role of 
(3 is played by J(0), and the density profile is given in terms of the magnetization 
variable m = 2p — 1. 

We are left with showing that there is no solution to (4.21) for h 7^ 0. In [9] (see 
also [11]) it is proven that the evolution equation 

d t u(r, t) = -u(r, t) + (f)- 1 (V * H; t)(r) - 1/2) + h) (4.22) 

has a unique solution in the class of travelling wave solutions if f'(p) — J(0)(p — 
1/2) — h = has three distinct solutions and this solution has non zero speed 
for h 7^ 0. In particular this implies that there is no solution with zero speed, i.e. 
stationary, in the case of three distinct roots. Due to (2.8), f'(p) — J(0)(p — 1/2) — h 
can have at most three roots. If it has only one root p±, then u must be constant, 
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since the asymptotic values of u must be equal to p\. We are then left with the 
case in which there is a root p\ and a double root pi- With no loss of generality let 
us assume that P2 > p\ so that lim 2 ^ 00 u(z) = p2- Due to (2.8) we have that 

f'(u) - J(0)(u - 1/2) — h = J*(u — 1/2) - J(0)(u - 1/2) > (4.23) 

whenever u(z) G [u, pz], for some u < p2, and the inequality is strict for u(z) G 
[u, P2). If zq = inf{z : u(z) < P2} < 00, then a contraddiction is easily established 
by (4.23) with z = zq. If z = 00, since u is monotonic, (4.23) is valid for all 
z sufficiently large. However this implies that u is unbounded, which establishes 
the desired contraddiction and this completes the proof of uniqueness. We remark 
that uniqueness holds with no need of (2.8) if we restrict ourselves to establish 
uniqueness in the class of functions which have asymptotic at z = ±00 which are 
symmetric with respect to 1/2, since this imposes h = in (4.21). □ 

We just remark that the regularity of the istantonic solution can be easily read 
out from (4.21) in terms of the regularity of J and /. 

5. Formal asymptotic expansions 

The computations regarding Statement 1 and Statement 2 follow closely the 
scheme of R. L. Pego's work [29]. We start by expanding the density profile in the 
bulk (outer expansion) and then we will do the same close to the interface (inner 
expansion): we will finally match the values at the boundaries. 

We recall that (2.2), (2.3) and (2.7) are assumed. With abuse of notation, J(R) 
(R G M + ) will mean J(r), for \r\ = R. Until further notice, r = r q , any value of q 
(see (3.1) and (3.2)). 

The outer expansion: we write 

p e (r, r) = p (r, t) + ept(r, r) + e 2 p 2 (r, r) + . . . (5.1) 

a e (r, r) = a(p e (r, r)) = a (r, t) + eai(r, t) + e 2 a 2 (r, r) + . . . (5.2) 
p e (r, t) = n (r, t) + ep x (r, r) + e 2 p 2 (r, r) + . . . (5.3) 
By Taylor expansion, using (3.3), one obtains directly 

(T (r,r) = a(p (r,r)) (5.4) 
0-1 (r,r) = pia'(po) (5.5) 
o~2 = P20- (po) + Pi — - — (5.6) 

Po = f'(po) ~ J(0)po (5.7) 
Pi = f" (po)pi - i(0)pi (5.8) 

P2 = f'(po)P2 + f"§ - J(0)p 2 - ^Ap (5.9) 

in which J = (1/d) J J(r)r 2 dr and we recall that J(0) = J J(r)dr 
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The inner expansion: We want to expand the density profile p e in a small neigh- 
borhood of the interface T T , r > 0. One possible choice of the neighborhood is 

K = {r G R d : |0(r,r T )| < e a } (5.10) 

where a G (0, 1) (recall that <fi denotes the signed distance), let us set z = <f>(r, T T )/e. 
If e is sufficiently small (depending on the minimum radius of curvature of the 
interface) then for each r G Af e there is a unique x = x(r) G T r such that r = 
x + v{x)ez. Notice that 

i/(i(r)) = V^(r,r T ) (5.11) 
for any r G Af e . Moreover we define 

V(x) = V(x(r)) = |^(r, r) (5.12) 

Note that dfi/dr is independent of r G A/" e , as long as x(r) is held fixed. For r£jV e 
will look for expansions of the form 

p e (r, r) = p (z,r,T) + epi (z, r, r) + e 2 p 2 (^, t) + . . . (5.13) 

cr e (r, r) = a (z, r, r) + e5"i(z, r, r) + e 2 o- 2 (^, r,r) + . . . (5.14) 

/i e (r, r) = /z (z, r, r) + e/ii(z, r, r) + t 2 jx 2 {z, r,r) + . . . (5.15) 

in which any quantity w(^, r, r) is defined for every r in Af e and iu(z, r + cv, r) = 
10(2:, r, r), f = z/(x(r)), for all c G R such that r+cu G A/" e , so that v- J \/w(z, r, r) = 0. 
The following formulas will be useful 

Vp(r, r) = e _1 z/(r, r)d z p(z, r, r) + V r /5(z, r, r) (5.16) 

<9 T p(r, r) = 9 T p(^, r, r) + e _1 l/(r, r)d z p(z, r, r) (5.17) 
The inner expansion for a is identical to the outer one 

a = o-(po) (5.18) 

cri = a'(p )p 1 (5.19) 

a2=a // (po)p 2 + a / (po)y (5.20) 

More complicated is the computation of ftf. for the moment we give it only to first 
order 

Ao = /'(p )- J*Po (5.21) 
in which * stands for the convolution in the z variable, that is 

(J ★ w) (z, r, r) = / J(z — z')w(z\ x, r)dz' 
Jr 

and we recall that J is defined in (3.5). 
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The relation between T T and pi is given by the following front centering condition 
[29]: 

U'(z)[p(z,r,r) -U(z)]dz = (5.22) 



where p = £V e l pi and the coefficient of each power of e in (5.22) is required to 
vanish. 



Matching conditions: We will match the outer and the inner expansion of the 
chemical potential order by order by requiring formally that 

(p-o + €Hi + e 2 p 2 + ■ ■ ■ ) (r + ezv(r, r) , r) (p, + ep x + + ...)(z,r,t) (5.23) 

with z = e 1_a , a G (0, 1). By expanding the left-hand side of (5.23) we obtain (see 
[29]) 

^o( r > r ) = lirn fio(z,r,r) (5.24) 

z—>±oo 

(pf + zv ■ Wpq)(t, t) = pi(z, r, t) + o(l) as z — > ±oo (5.25) 
+ 2z/ • V/U^ + 2^ 2 ( z/ ' ^) 2 A t o = ^ ( r i r ) = ^(-z, ic, t) + o(l) as z — > ±oo (5.26) 
in which jtt^(r, r) = lim t ^ ± pi(r + vt, r), and analogous formula for the derivatives. 

Verification of Statement 1: the case q = 2. The starting point is equation (3.2) 
and we set r = r^. By using (5.11), (5.16) and (5.17), we obtain that the highest 
order arising in the inner expansion of (3.2) (0(e~ 2 )) yields the equation 

d z (a d z p ) = 0. (5.27) 

Using now the explicit expressions for a (5.18) and p (5.21) we see that the 
solution p of (5.27) coincides with the stationary solution of (3.4), that is 

p (z,r,r) = U(z) (5.28) 

Moreover by using (3.6) we obtain 

Ao(^,r,r) = ^ (5.29) 
The leading order of (3.2) in the outer expansion (0(1)) is 

d T po = V (a Vpo) (5.30) 

which by the assumptions on (5. 4), (5. 7) (see (2.7)) and on the initial condition is 
a nondegenerate parabolic PDE. (5.30) has to be supplemented with the boundary 
conditions coming from the matching rule (5.24) and the asymptotic value of po, 
given by (3.7). This completes the derivation of (3.11) and we are left with the 
identification of the speed of the interface. To do this we have to look at the next 
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order in the inner expansion (0(e 1 )). From (5.12), (5.16), (5.17) and the fact that 
jlo is constant (5.29) we obtain 

Vd z p = d z (dodofii) (5.31) 

Recall that po = U and integrate in z E K. to obtain 

« (8 32) 

in which = iu(+oo, r, r) — u;(— oo, r, r) and recall that cr(p + ) = cr(p~) = a^. 

Finally the matching condition (5.25) implies 

[dzM + -Z = [vVno£ (5.33) 

and the notation [•] 1 is introduced in Statement 1. Formula (5.32) together with 
(5.33) implies (3.12). This completes the verification of Statement 1. □ 



Verification of Statement 2: The case q = 3. 

In this section r = T3. By assumption p (r, 0) = p ± in and the the highest 
order in the outer expansion of (3.2) implies that po(r,r) = p ± in ttf also for 
r > 0. The lowest order in the inner expansion is still 0(e -2 ) and it yields once 
again (5.27), so that po = U, and (5.29). The next order in the outer expansion 
(0(e)) gives 

d T Po = V (a V^i + arVpo) (5.34) 
but since d T p = 0, Vpo = and a Q (p ± ) = a ± 7^ we obtain 

A^i = (5.35) 

Since d z po = (5.29), the 0(e _1 ) in the inner expansion gives 

d z {a d z p 1 } = (5.36) 

so that d z pi = 0, since ao 7^ and by (5.25) pi(z) approaches (as z approaches 
±00) the value of v • Vpo near the interface and Vpo = 0. Hence 

~p Jl (z,r,r) = C(r,r). (5.37) 

We will verify below that 

fa(z,r,r) = pi(z, r, r)f"(p (z, r, r)) - J * p x - K(x(r)) J J(z - z')p (z' ,r,r)dz' 

(5.38) 

and putting together (5.38) and (5.37) we obtain 

~ Pl (z, r, r)f"(U(z)) - (J* ~ Pl )(z, r, r) - K(x(r)) [ (z - z')J(z - z')dz' = C(r, r) 

Jr 

(5.39) 
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Multiply both sides of (5.39) by U'(z) and integrate in the variable z over all M. 
Since J pi[d z f'(U)/ — J*U']dz = 0, the expression between square brackets being 
zero because it is simply <9 2 //o, we obtain 

-K(r) [ (z- z')J{z- z')U{z')U'{z)dzdz' = C{r,r)[U}±™ (5.40) 

JRxR 

In the left-hand side of (5.40) we recognize the surface tension of the model (3.13) 
and so 

fn = C = (5.41) 

and in particular fli is independent of time. The velocity of the interface is once 
again retrieved by looking at the next order in the inner expansion (0(1)) 

VU' = d z (a d z p 2 ) (5.42) 

so that by integration 



a± [<9^ 2 ] + °° 

v = 

and by (5.26) we finally obtain 



V = " ^~7- J -°° (5.43) 



V 3 = V = 1 , J ~ 5.44 

[u] + T 



-oo 



that is (3.14). On the other hand (5.25) and (5.41) give us the boundary conditions 
for fix, which then solves A^i = with 

/Xl ( r ) = _£(rg (5.45) 

when r G T r , and the boundary value problem (3.15) is derived. We are left with 
showing that the first order correction to the chemical potential in the bulk is 
given by (5.38). The expression in the right-hand side of (5.38) consists of three 
terms: the first one is simply /"(po)pi- The other two originate from the 0(e) 
contributions to the integral term 

Jje(r- r') [po (0(r', T T )/e, r', r) + e~ Pl (0(r', T T )/e, r', r) + 0(e 2 )] (5.46) 

The 0(1) contribution to (5.46) is simply J * p and if we subtract it from (5.46) 
and we use po = U we obtain 



J J e (r -r')U(z')dr' - (j : *U) (z) 



+ eJ*pi + 0(e 2 ) (5.47) 



in which z (respectively z') is the rescaled distance of r (respectively r') from the 
interface, as before, i.e. z = <f>(r,T T )/e and z' = <p(r' ,T T )/e. We are then left 
to show that the term between square brackets is of order e and that its value 
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divided by e coincides to leading order with the third term in the right-hand side 
of (5.38). By choosing an appropriate frame of reference, we can assume that 
Ti = for i = 1, . . . , d — 1 and r d = ze and that in this coordinate system the 
upper half plane {r G W 1 : ra = 0} is tangent to T T . We will also orient the i-axis, 
i = 1, . . . , d — 1 along the directions of principal curvatures. The curvature in the 
direction i will be denoted by ki, so that K = Yli=i ki- It is moreover convenient 
to go back to the original coordinates, that is to x = re -1 = (0, . . . , 0, z) and 
x' = r'e -1 : z' = <j)(r',T T )/e is left unchanged, since 4>{x' , e _1 r r ) = (p(ex\T T )/e. 
We have 

Notice that since J is supported in £?i(0),the <i-dimensional ball of radius 1, cen- 
tered at 0, we are only interested in \x'\ = O(l) and z' = O(l). We will pa- 
rametrize e _1 r r as a function f e : M. d — > R around the point G M d_1 , that is 
{r G R d : |(n,... 7 r d - 1 )\}ne- 1 T T = {(qj e (q)) : q G R d ~\ \q\ < 1}. By the fact 
that T T is assumed to be smooth we have 

/ e (9) = -E^ + 0(e 2 ) (5-49) 
i=i 

for \q\ < 1. Our aim is to write z' explicitely as a function of x' so that it will make 
the integral in (5.48) more explicit. Since x' G e _1 A/" e we have 

x' = (q, fe(q)) + -fc==^' + 0(e 2 ) (5.50) 

v 1 + |v/ e r 

and, by using (5.49), for i = 1, . . . , d — 1 we obtain 

x' i = q i + ehqiz' + 0(e 2 ) (5.51) 



and 



which imply 



x' d = z'-eJ2^f + 0(e 2 ) (5.52) 



i=i 



d-l , ,2 

rvjX a 



z' = x' d + eY,-^ + 0^) (5.53) 

i=i 

Let us now go back to (5.47). The term in the square brackets is equal to 

J J(x - x') [U(z') - U{x' d )] dx' = e J J(x - x') U'{x' d ) dx' + 0(e 2 ) 

(5.54) 

in which we used the definition of J (3.5) and (5.53). We are then left with evalu- 
ating 

/ J(x ~ x') ^ j u> {x > d)dx > (5.55) 
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in which x = (0,... ,0,z). Set q = (x[, . . . ,a^-i) G = z' and 

vV + (z -z') 2 . Then the expression in (5.55) is equal to 



lH k iJ J(R)U'(z') q fdqdz' 



'd-\ 



= 1(11^) J J(R)U'(z') q ldqdz' 
= -f / (^) ^)<7i 2 cW (5.56) 

= -f / gi J(fl)(^-*M^)?idgd*' 

= y/ wc*'-*)^*')^' 

in which the first equality follows because J is rotation invariant (isotropy), the 
second is integration by parts with respect to z' and the fact that K = h, the 
third is the chain rule of differentiation, the fourth is integration by parts with 
respect to q± and the last one follows from the definition (3.5). 
This ends the verification of Statement 2. □ 

Appendix A: the surface tension. 

This Appendix generalizes [4], which deals with the case of / given in (2.14). We 
are in the same framework as Proposition 4.2. The surface tension can be defined 
as the difference between the free energy of an equilibrium state with an interface 
and a homogeneous one [4,33], i.e. from (2.1) 

^ r+L p+L j-+M 

S = L—*oo (2Lp^ ^ J L dXl ---J L dxd - 1 J M dXd - /c( ^ +) ] 

(6.1) 

in which p*(x) = U(ed • x), where = (0, . . . , 0, 1) G M d , and for any measurable 
function p : R d -> [0, 1] 

g(p)(x) = f(p(z)) - \[p{x) - 1/2](J * [p - l/2])(x) (6.2) 

Observe that f c (p + ) = fc(p~) = 9(p + ) = 9(p~)- (6-1) is well defined because of 
Proposition 4.2 and the properties of /. Clearly (6.1) reduces to 

S = [ [g(p*)(0,...,0,z)-f c (p+)]dz. (6.3) 
Jr 

We are going to show that S, as defined in (6.1), coincides with (3.13). 
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Observe that for p e C 1 



" (P) - f\P)^ -\&'\P- V2] " i[P " 1/2] J * £ (M) 



but, by (3.6), /'(p*) = J * (p* - 1/2). Hence 

= ^ (u'(x d )(j*[U - l/2])(x d ) - \U(x d ) - l/2](.7* t/')M) 
Integrating (6.3) by parts and substituting (6.5) in the integrand we finally obtain 

s = ~\ I Xd-p-(p*)dx d 
2 Jr ox d 

= \ I [( z '~ z)U'{z)J{z - z')[U{z') - 1/2]<W 
2 Jr Jr 

= lf f \z' -z)U'(z)J(z-z')U(z')dzdz' (6.6) 
2 Jr Jr 

where we used the fact that there is an m' > such that \im z ^± 00 e m \ z \U'(z) = 0, 
which follows from (3.6) and Proposition 4.2, point (3). Since (6.6) is (3.13), we 
are done. 

ACKOWLEDGEMENTS 

It is a pleasure to acknoweledge several fruitful discussions with P. Butta, E. 
Orlandi, R.V. Kohn, E. Presutti, L. Triolo and H. Spohn. G.G. would like to 
thank B. Brighi and M. Chipot for several discussions and considerable help on the 
material in Section 4. This work was partially supported by the Swiss National 
Science Foundation, Project 20-41'925.94 (G.G.) and by NSF-DMR 92-13424 4- 
20946 and AFOSR 0159 4-26435 (J.L.L.). 

References 

[1] N. Alikakos, P. W. Bates, and X. Chen, Convergence of the Cahn-Hilliard equation to the 

Hele-Shaw model, Arch. Rat. Mech. Anal. 128 (1994), 165-205. 
[2] L. Bonaventura, Interface dynamics in an interacting spin system, Nonlinear Analysis TMA 

25(8) (1995), 799-815. 
[3] L. Bronsard and R.V. Kohn, Preprint (1994). 

[4] P. Butta, On the validity of a Einstein relation in models of interface dynamics, J. Stat. 
Phys. 72 (1993), nos. 5/6. 

[5] G. Caginalp, The dynamics of a conserved phase field system: Stefan like, Hele-Shaw and 
Ginzburg-Landau dynamics, IMA, J. Apll. Math. 44 (1990), 77-94. 

[6] J.W. Cahn, CM. Elliott and A. Novick-Cohen, The Cahn-Hilliard equation: motion by 
minus the Laplacial of the mean curvature, Euro. J. Appl. Math. 7 (1996), 287-301. 

[7] J. W. Cahn and J. E. Hilliard, Free energy of a nonuniform system. I. interfacial free energy, 
J. Chem. Phys. 28 (1958), 258-267. 

[8] X. Chen, Hele-Shaw problem and area-preserving curve shorting motion, Arch. Rat. Mech. 
Anal. 123 (1993), 117-151. 

[9] X. Chen, Existence, uniqueness, and asymptotic stability of traveling waves in nonlocal evo- 
lution equations, Preprint (1995). 

22 



R. Dautray, J.-L. Lions, Analyse mathematique et calcul numerique pour les sciences et les 
techniques, tome 3, Masson, Paris (1985). 

A. De Masi, T. Gobron and E. Presutti, Travelling fronts in non-local evolution equations, 
Arch. Rat. Mech. Anal. 132 (1995), 143-205. 

A. De Masi, E. Orlandi, E. Presutti, and L. Triolo, Glauber evolution with Kac potentials: I. 
mesoscopic and macroscopic limits, interface dynamics, Nonlinearity 7 (1994), 1-67. 
A. De Masi, E. Orlandi, E. Presutti, and L. Triolo, Stability of the interface in a model of 
phase separation, Proc. Royal Soc. Edinburgh 124A (1994), 1013-1022. 

A. De Masi, E. Orlandi, E. Presutti, and L. Triolo, Uniqueness and global stability of the 
instanton in non local evolution equations, Rendiconti di Matematica, Serie VII 14 (1994), 
693-723. 

A. De Masi and E. Presutti, Mathematical methods for hydrodynamic limits, Lecture notes 
in mathematics. Springer- Verlag, Berlin 1501 (1991). 

P. De Mottoni and R. Dal Passo, The heat equation with a nonlocal density dependent ad- 
vection term, unpublished manuscript (1991). 

L. C. Evans, H. M. Soner, and P. E. Souganidis, Phase transitions and generalized motion 
by mean curvature, Comm. Pure Appl. Math. 45 (1992), 1097-1123. 

G. Giacomin and J.L. Lebowitz, Exact macroscopic description of phase segregation in model 
alloys with long range interactions, Phys. Rev. Lett. 76(7) (1996), 1094-1097. 

G. Giacomin and J.L. Lebowitz, Phase segregation dynamics in particle systems with long 
range interactions I: Macroscopic limits, J. Stat. Phys. to appear (1996). 
D. Gilbarg, N.S. Trudinger, Elliptic partial differential equations of second order (second 
edition), Springer- Verlag (1983). 

M. Kac, G. Uhlenbeck, P.C. Hemmer, On the Van der Waals theory of vapor-liquid equilib- 
rium. I. Discussion of a one dimensional model, J. Math. Phys. 4 (1963), 216-228. 
M.A. Katsoulakis, P.E. Souganidis, Interacting particle systems and generalized evolution of 
fronts, Arch. Rat. Mech. Anal. 127 (1994), 133-157. 

O.A. Ladyzenskaja, V.A. Solonnikov and N.N. Uralceva,, Linear and quasilear equations of 
parabolic type, AMS, Translations of Mathematical Monographs, Providence, Rhode Island 
(1968). 

J. S. Langer., An Introduction to the Kinetics of first-order phase transitions, Solids far from 
equilibrium. C. Godreche ed., Cambridge Univ. Press (1991). 

J. L. Lebowitz, GHS and other inequalities, Comm. Math. Phys. 35 (1974), 87—92. 
J. L. Lebowitz and O. Penrose, Rigorous treatment of the Van der Walls-Maxwell theory of 
the liquid-vapour transition, J. Math. Phys. 7 (1966), 98. 

J.-L. Lions, Equations differentielles operationelles et problemes aux limites, Springer- Verlag 
(1961). 

National Reserarch Council., Mathematical research in materials science: opportunities and 
perspectives, Board on Mathematical Sciences. National Academy Press, Washington, D.C. 
(1993). 

R. L. Pego, Front migration in the nonlinear Cahn-Hilliard equation, Proc. R. Soc. Lond. A 
422 (1989), 261-278. 

O. Penrose and J. L. Lebowitz, Rigorous treatment of metastable states in the Van der 
Walls-Maxwell theory, J. Stat. Phys. 3 (1971), 211-236. 

B. Simon, The statistical mechanics of lattice gases (I), Princeton University Press (1993). 

H. Spohn., Large scale dynamics of interacting particles, Texts and Monographs in Physics. 
Springer- Verlag, Berlin (1991). 

[33] H. Spohn, Interface dynamics in particle models, J. Stat. Phys. 71(5/6) (1993), 1081-1131. 



GlAMBATTISTA GlACOMIN 

institut fur angewandte mathematik 
Universitat Zurich-Irchel 

wlnterthurerst. 190, ch-8057 zurich switzerland 
E-mail address: gbg@amath.unizh.ch 



Joel L. Lebowitz 
Department of Mathematics and Physics 

23 



Hill Center, Rutgers University 
New Brunswick, N.J. 08903, U.S.A. 

E-mail address: Lebowitz@math.rutgers.edu 



24 



